o 
o 

(N 

< 



o 
o 

u 

Oh' 



I 

O 

o 



> 

O 

o 
o 

-I— > 



I 

o 



X 



Survival of superconducting correlations across the 2D superconductor-insulator 

transition: A finite frequency study 
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The complex AC conductivity of thin highly disordered InOi^ films was studied as a function of 
magnetic field through the nominal 2D superconductor-insulator transition. We have resolved a 
significant finite frequency superfluid stiffness well into the insulating regime, giving the first direct 
evidence for superconducting correlations in the insulating state of an amorphous film. A phase 
diagram is established that includes the superconducting state, a transition to a 'Bose' insulator and 
an eventual crossover to a 'Fermi' insulating state at high fields. We speculate on the consequences 
of these observations and their impact on the understanding of the insulating state and purported 
intervening low temperature metal. 
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A quantum phase transition (QPT) is a zero temper- 
ature change of state as a function of some non-thermal 
parameter (pressure, doping, magnetic field, etc.)i^. The 
2D superconductor-insulator quantum phase transition is 
a particularly beautiful and illustrative one, reflecting a 
transition between the two eigenstates at the extremes 
of a superconductor's fundamental uncertainty relation 
between phase and particle number (A^An > 1). One 
of the key questions here is the nature of the zero tem- 
perature destruction of the superconducting state. Does 
it proceed in a mean-field fashion by destruction of the 
amplitude A of the superconducting order parameter 
^p = Ae'^ or is it dominated by fluctuations of the su- 
perconducting phase 97 The manner in which the su- 
perconductivity is destroyed is an issue that has direct 
relevance to many other important problems including 
that of high-temperature superconductivitj*^. 

A related matter is the nature of the insulating state. 
There have been a number of proposals for an insulator 
with strong superconducting correlations. For instance, 
Fisher and co-workers^ postulated a dual description - 
the so-called 'dirty boson' model - of the SIT in which 
the superconducting state reflects the condensation of 
Cooper pairs and localization of vortices, while the in- 
sulating state is characterized by condensed vortices and 
localized Cooper pairs. More recently, it has been pro- 
posed that the transition is an inhomogeneous one&, and 
where global superconductivity may obtain by percola- 
tion of locally superconducting clusterai. To date, ev- 
idence for superconducting signatures in the insulating 
state may have been found in signs of vortex activation 
in insulating fllma^ as well as in a crossover of the Hall co- 
efficient Rx-jj^ at a fleld higher than the superconductor- 
insulator critical field Us it- Although these results have 
been suggestive there has been a notable lack of direct 
evidence for superconducting correlations in the insula- 
tor. Moreover a number of outlying issues remain. For 



instance, the Cooper pair gap appears to close on the 
approach to the SIT in amorphous films, evidence seem- 
ingly incompatible with the existence of localized Cooper 
pairsifliii. Still others have challenged the existence of a 
direct transition between superconducting and insulating 
states altogether and instead postulate the existence of 
an intervening metaU^. 

In this Letter we present the results of a comprehensive 
study of the complex AC conductivity through the nom- 
inal 2D SIT in InO^ thin films using microwave cavities. 
Here, we have resolved a significant finite frequency su- 
perfiuid stiffness in the insulating state which persists up 
to resistances over 10^ Ohms/D and fields up to 3x the 
critical field, giving the first direct evidence for a state 
with localized Cooper pairs. This establishes a phase di- 
agram with distinct regions dominated by superconduct- 
ing, 'Bose' insulating, and 'Fermi' insulating effects. 

Samples were 200 A-thick 3mm-diameter highly- 
disordered amorphous indium oxide (a:InO) thin films 
prepared by e-gun evaporating high purity (99.999%) 
In2 03 on clean 19mm-diameter sapphire discs in high 
vacuumiii^. Post-deposition, samples were room- 
temperature annealed for approximately one week in am- 
bient air and then, except for short periods, were kept 
well below LN2 temperatures for the duration of the 
measurements (approximately 2 months). TEM mea- 
surements show diffuse rings with no diffraction spots, 
suggesting the amorphousness of the films with no crys- 
talline inclusions. AFM images show continuous films 
with no voids or cracks and are completely featureless 
down to a scale of a few nm (the resolution of the AFM) . 
Experiments were performed in a novel cryomagnetic 
resonant microwave cavity system. The cavity diame- 
ter was chosen to optimize performance in the 22 GHz 
{hui/kB = 1.06 K) TEOll mode, but measurements were 
possible at a number of discrete frequencies from 9 to 
110 GHz. Cavities were operated at a power low enough 
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FIG. 1: (a) DC zero field sheet resistance vs. temperature, 
together with a fit to the AL form with mean-field transi- 
tion temperature Tco = 2.16K. (b) DC Sheet resistance vs. 
temperature for field values H=0, 0.05, 0.1, 0.25, 0.5, 1, 1.5, 
2.5, 3, 3.5, 4, 4.5, 5, 5.5, 6, 6.5, 7, 7.5, 8 Tesla. (c) Sheet 
resistance vs. magnetic field for temperatures shown. The 
crossing point determines a critical field of 3.68 Tesla. 



to prevent sample heating. Various well-known relations 
were used to relate the resonances' frequency shift Alu 
and change in quality factor AQ upon sample introduc- 
tion to the complex conductivitjii^*iSiiLi&. The conver- 
sions from AQ and Alu to conductivity were made by 
normalizing the AC resistance of the cavity data to DC 
data at temperatures well above the occurrence of su- 
perconductivity and insisting that the superfluid stiff- 
ness (defined below), was frequency independent at zero 
field as T — > 0. DC resistance was measured on co- 
deposited samples in a two-probe configuration by low 
frequency AC lock-in techniques using excitation cur- 
rents of lOpA — lOnA. The probe's lead resistances, 
which have a negligible temperature dependence in the 
displayed temperature range, were well characterized and 
have been subtracted from the displayed data. 

The zero field DC resistance curve shown in Fig. la 
is fairly typical for a highly disordered superconducting 
thin film showing an approximately log increase of R 



with decreasing temperature, before the occurrence of a 
broad superconducting transition. As shown in Fig. lb, 
when the field is increased the resistance curves cross over 
from a superconducting positive dR/dT to an insulating 
negative dR/dT behavior similar to previous field-tuned 
studiesi^iSaiSi. When plotted asRvsH (Fig. Ic) the ex- 
perimental data shows a low-temperature iso-resistance 
crossing point of 3.68 Tesla, which can be identified as 
the critical field of the SIT, Hsit- The DC data can 
be shown to be consistent with previous studiesi^iSSiSi 
that found scaling as a function of the reduced variable 
\H — Hsit\/T^^^'^ where zv is consistent with the expo- 
nent for 2D classical percolation |. 

As the temperature is lowered, distinct regions typi- 
fied by amplitude and then phase fluctuations of the su- 
perconducting order parameter ip = Ae are expected. 
We can quantify the amplitude fluctuations at zero field 
by fitting the DC resistance curves to the Aslamazov- 
LarkinS^ (AL) form for gaussian fluctuations. In two di- 
mensions, the extracted temperature scale Tco = 2.16K 
does not signify the occurrence of a phase transition, but 
instead represents the temperature scale below which the 
Cooper pair amplitude becomes well defined and as such 
the AL form is not a good fit near TcO itself. 

We now turn our attention to the AC conductivity. 
At low temperatures, the imaginary conductivity for a 
long-range ordered superconductor is expected to have 
the form (72 = -^^^ where N is the superfluid density 
and e and m are the electronic charge and mass respec- 
tively. For a fluctuating superconductor one can define 

o'2 — — where an additional frequency dependence is 
captured by a generalized frequency dependent superfluid 
density. The superfluid density is directly proportional 
to the superfluid stiffness which is the energy scale for 
inducing slips in the superconducting phase. In Fig. 2a 
we display the H = generalized frequency dependent 
superfluid stiffness, Tg (in degrees Kelvin) extracted via 

the relation (72 = o-q ^f^" , where ctq = ^ is the quan- 
tum of conductance for Cooper pairs divided by the film 
thickness. The superfluid stiffness curves all cross the 
predictive line {Tg — ATktb) for the Kosterlitz-Thouless- 
Berezinskii (KTB) transition in a range near 1.6 K. We 
associate this temperature as the KTB scale below which 
the superconductor becomes robust against vortex phase 
fluctuations. A somewhat similar phenomenology has 
been seen in quasi 2D high temperature superconduc- 
tors at THz frequencies^, although here the superfluid 
stiffness acquires a frequency dependence below the KTB 
scale. The Tktb of 1.6 K and T^o of 2.16 K give the 
temperatures at zero field for phase coherence and a well 
defined superconducting amplitude, respectively. 

Displayed in Fig. 3 is the finite-frequency complex re- 
sponse at i7 ^ for the 22 GHz data. Similar plots can 
be made at other frequencies. Fig. 3a shows the AC 
resistance measured over a range of temperatures and 
magnetic fields. We identify the region below the critical 
contour of i? = 2295ri (in bold) as the domain over which 
the underlying T = superconducting phase has influ- 
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FIG. 2: Superfluid stiffness, Tg (oc u>a2), extracted as de- 
scribed in the text, (a) Temperature dependence at H=0. 
Also shown is the theoretical Tktb which intersects Tg around 
T = l.QK. (b) 0.5 K magnetic field dependence at select low 
frequencies. 



ence. In Fig. 3b is shown the superfluid stiffness Tq at 22 
GHz, using the same conversion as above, and plotted as 
a function of temperature and magnetic field. The use of 
relatively high probing frequencies allows us to resolve su- 
perconducting fluctuations into the insulating part of the 
phase diagram. An additional advantage of our high fre- 
quencies is that they are at least two orders of magnitude 
higher than a generous estimate for the vortex depinning 
frequencji2& found in thin films of conventional supercon- 
ductors. The vortex contribution to the optical response 
is therefore expected to be purely dissipative, leaving the 
superfluid as the only principal contributor to cr2^- The 
finite-frequency superfluid stiffness falls quickly with in- 
creasing fleld, but remains flnite above Hsit, well into 
the insulating regime to flelds almost 3 times the crit- 
ical field H5/T. To the best of our knowledge this is 
the first direct measure of superconducting correlations 
on the insulating side of the 2D superconductor- insulator 
transition in an amorphous film. 

Our observation of a flnite frequency superfluid stiff- 
ness Sit H > HsiT is not inconsistent with an insulating 
T=0 groundstate. As alluded to above, our experiments 
are sensitive to superfluid fluctuations because we probe 
the system on short time scales via an experimental fre- 
quency uJExp that is presumably high compared to an in- 
trinsic order parameter fluctuation rate ujqc close to the 
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FIG. 3: (a) AC Sheet resistance (« Re{-^)) and (b) Super- 
fiuid stiffness (oc uja2) at 22 GHz. The critical field, Hsit, is 
shown as a white dot at H = 3.68T. Yellow indicates maxi- 



transition. Above Hsit, an experimental probe in the 
limit ujExp ^ '-^QC will not detect superfluid fluctuations 
as a system with Tg{uj —>■ 0) ^ can support superfluid 
flow, behavior obviously incompatible with the notion 
of an insulator. On general grounds an insulator with 
(Ti = at zero frequency, must have a (72 that is negative 
a.t uj = Q'^ by Kramers-Kronig considerations, meaning 
that such a system can not appear superconducting at 
low frequencies. It is only by using a relatively high uj 
probe we detect these fluctuations. Returning to Fig. 
2b we see the dramatic drop in superfluid stiffness with 
increasing field and the strong frequency dependence of 
Tg which refiects that the Cooper pairs are correlated 
on time scales that can be resolved at the probing fre- 
quency, so that while long range order does not exist at 
finite T, short-range correlations do. We also note that at 
low temperatures and well into the insulating side of the 
phase diagram, that although the magnitude of the su- 
perfluid stiffness is of the same order as it above Tco, the 
signal becomes temperature independent at low temper- 
ature showing the intrinsic quantum mechanical nature 
of these fluctuations. 

We are able to extract a phase diagram that estab- 
lishes the existence of superconducting correlations well 
into the insulating state. In Fig. 4 contours are plotted 
which denote the region above which our superfluid stiff- 
ness becomes almost indistinguishable from the normal 
state noise level (set at 1 % of the T — > 0, iJ = super- 
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FIG. 4: 2D Field Tuned Superconductor-Insulator "Phase Di- 
agram" . Contours with markers are the detection limit for Te 
at specified frequencies. Black dots on the temperature axis 
denote Tamv (= ^co) and Tphas^ (= Tktb), temperatures 
signifying the onset of amplitude and phase fluctuations, re- 
spectively. HsiT appears as a black dot and the contour of 
critical resistance from Fig. 3a appears as a solid black line. 
Open symbols represent data obtained from a small linear 
extrapolation beyond our maximum field of 8 Tesla. 



fluid stiffness), thereby giving a measure of the extent of 
superconducting correlations into the insulating regime. 
The noise is greatest for frequency contours away from 
our cavity's optimal operating frequency of 22 GHz. It is 
evident that the higher frequency probe at 22 GHz allows 
one to examine the fluctuations of the order parameter 
persisting at fields higher than Hsit and temperatures 
higher than Tco, at shorter length and time scales than 
the 9 or 11 GHz probes. In the high-frequency limit one 
would expect that the detection limit would eventually 
extrapolate to a field where the Cooper pairs are com- 
pletely depaired. For 22 GHz we observe that the data 
extrapolates to a value of H=9.35 T, which is close to the 
value of the pair breaking field, ifc2, found in similar films 
of InOafi^iSL2& giving evidence that 22 GHz is near the 
high-frequency limit. It is impossible to say with our cur- 
rent results whether this pair breaking scale should be as- 
sociated with orbital or spin effects. We do note that this 
field is well above the naive estimate [Hp = 1.86 x Tc) of 
the Pauli pair breaking field, but that disorder is known 
to lead to enhanced Hc2 over the expectation from its 
H=0 valuers. Our results and interpretation are not con- 
sistent with the inference of RefiSi that Cooper pairs still 
had integrity at fields many times the pair breaking scale. 



It would interesting to verify whether the crossover ob- 
served in Rxy by Paalanen et al^ - which has also been 
interpreted as a 'Bose' glass/'Fermi' glass crossover - is 
found in our films at the extrapolated pair breaking field 
scale. 

We emphasize that we believe the only true phase tran- 
sition in the system (neglecting for a moment the pos- 
sibility of an intervening metal phase) is at Hsit and 
that the contours in Fig. 4 give a frequency dependent 
crossover. We observe superconducting correlations not 
just asymptotically close to the SIT, but an extended re- 
gion above Hsit- This establishes for the first time a 
model-free, unambiguous picture of an insulating state 
that is dominated by superconducting correlations - a 
'Bose' insulator. Our phase diagram is characterized by 
a superconducting groundstate, a phase transition to a 
'Bose' insulating state at Hsit and a crossover to 'Fermi' 
insulator near a depairing field "Hc2" ■ 

This work raises certain questions about the expected 
electrodynamic response of the various proposed phases 
of matter. It might be expected, for instance, that the 
functional dependence of the DC R{T) is different for 
'Fermi' vs 'Bose' insulators, but these differences may be 
subtle and difficult to distinguish experimentally. In the 
present case, we have a heuristic argument about the de- 
tection of finite frequency superfluid stiffness as being in- 
dicative of the 'Bose' insulator, but detailed calculations 
are lacking. Moreover, it has been argued that an anoma- 
lous metallic state intervenes between the superconduct- 
ing and insulating states at low temperaturoiS. At finite 
temperature the difference in the dissipative response be- 
tween a fluctuating superconductor and this metal could 
be subtle, but the difference in the reactive response may 
be dramatic. Again, detailed calculations are lacking. In 
either case we have shown that the underlying supercon- 
ductivity will manifestly have to be taken into account 
to describe the insulating and anomalous metallic states. 
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